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The low energy band structure of the FeAs based superconductors is fitted by a tight binding 
model with two Fe ions per unit cell and two degenerate orbitals per Fe ion. Based on this, su- 
perconductivity with extended s-wave pairing symmetry of the form cos k x + cos k y is examined. 
The local density of states near an impurity is also investigated by using T-matrix approach. For 
the nonmagnetic scattering potential, we found that there exist two major resonances inside the 
gap. The height of the resonance peaks depends on the strength of the impurity potential. These 
in-gap resonances are originated in the Andreev's bound states due to the quasiparticle scattering 
between the hole Fermi surfaces around T point with positive order parameter and the electron 
Fermi surfaces around M point with negative order parameter. 
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The recent discovery of a new family of superconduc- 
tors, i.e. the FeAs based superconductors [1-6], has at- 
tracted much attention in the condensed matter com- 
munity. It has been reported that the superconducting 
transition temperature T c can be obtained as high as 55K 
[2] . The undoped iron arsenides have a spin density wave 
order below 150K [4]. When holes or electrons are doped, 
the iron arsenides become superconducting. 

Similar to the cuprate superconductors, the FeAs 
based superconductors also have a layer structure. It 
has been accepted that superconductivity comes from 
the Cooper pairs in the Fe-Fe plane. However, in the 
FeAs based superconductors, each unit cell contains two 
Fe ions and two As ions. The four As ions around each Fe 
ion do not locate in the Fe-Fe plane and have a two-fold 
rotation symmetry and two reflection symmetries (see 
Fig. 1). Due to different arrays of As ions around Fe 
ions, the Fe-Fe plane can be divided into two sublattices 
A and B. We note that the diagonal directions of the Fe- 
Fe plane have the translational symmetry with the period 
a. In this coordinate system, the momentum is a good 
quantum number. 

Angle resolved photoemission spectroscopy (ARPES) 
experiments have probed the electronic properties in the 
FeAs-based superconductors [7-15]. It is established that 
there are two hole Fermi surfaces around (0,0) and two 
electron Fermi surfaces around (7r,7r). These Fermi sur- 
face characteristics have been obtained by the LDA cal- 
culations [16-19]. Many theoretical models have been 
presented to reproduce the hole and electron pockets by 
employing the Fe d and As p orbitals and the hybridiza- 
tion among them [20-25]. However, there is no consensus 
on the superconducting gaps on the Fermi surfaces. In a 
series of ARPES, scanning tunneling microscopy (STM) 
experiments, and point-contact Andreev reflection spec- 
troscopy experiments, the order parameter has been in- 
terpreted to be nodeless [7,9-15,26], nodal [27-29], single 
gap [10,26,29-32], and multiply gaps [7,11-15]. 
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FIG. 1: (Color online) Schematic lattice structure of the FeAs 
layers with each unit cell (denoted by i and j) containing two 
Fe (A and B) and two As (A and B) ions. The As ions A 
and B are located just above and below the center of each 
face of the Fe square lattice, respectively. Here, ti is the 
nearest neighboring hopping between the same orbitals d xz 
or d yz . t-2 and are the next nearest neighboring hoppings 
between the same orbitals mediated by the As ions B and 
A, respectively. t4 is the next nearest neighboring hopping 
between the different orbitals. 



In this work, we start from two Fe ions per unit cell 
and two degenerate orbitals d xz and d yz per Fe ion and 
construct an effective four-band model, which exhibits 
the features of the Fermi surfaces in the FeAs-based su- 
perconductors. Based on the mean field theory for super- 
conductivity, we investigate the differential conductance 
and the impurity effect for the extended s-wave pairing 
symmetry [12,17,20], so that we can understand the elec- 
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FIG. 2: (Color online) The band structure of the four-band 
model with ti = 0.5, t 2 = 0.2, t 3 = -1.0, U = 0.02, and 
|U = -0.622 (eV), plotted along the path (0,0) -» (tt,0) -v 
(jr,jr)->(0,0). 



tronic properties in the FeAs-based superconductors. 

We assume that ti is the hopping between the same 
orbitals on the nearest neighboring Fe sites, t% and t% 
are the next nearest neighboring hoppings between the 
same orbitals mediated by As ions B and A, respectively, 
and ti is the hopping between the different orbitals on 
the next nearest neighboring Fe sites (see Fig. 1). It is 
expected that t± is small and has the same value in both 
translation symmetry directions. Therefore, the model 
Hamiltonian we propose can be written as 

Hq = — ^ {f 1 { c Aa.ii<T C Aa,i3<T + c Ba,ija C Bu,ij<j) 



+ [tl c 'Aa,ijij( C Ba,ija + CBa,i+lja + CBa,ij + la + CBa,i+lj + lcr) 
+ t 2{c\ alj(J CAa,i+l]a + c B a ,i]<? C Bu,ij + lcr) 
+ t 3(c\ aijcr C Aa ^j + la + C Baij(T CBa,i+ljo) 
+ti{ c Au,ija C Aa+l,i+l ] o + c\ aij(J CAa + l.i] + la 
+C~B a ,i]a C Ba+l,i+lj a + CBa,ijtr C Ba+l,ij+l<r) + h.C.]}, (1) 

where a is the spin index, i, j label the position of unit 
cell, and a = and 1 represent the degenerate orbitals 
d xz and d yz , respectively. Obviously, H possesses the 
same symmetry with the FeAs-based superconductors, 
which is key to understand the electronic properties of 
this new family of high temperature superconductors. 

By diagonalizing the above Hamiltonian in the momen- 
tum space, one obtains 
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where u(v) — 0, 1, e^.k = —2(^2 cos k x +t% cos k y ),es,k = 
—2(t 2 cos k y + t 3 cos k x ), e xy> k — — 2^4(cos k x +cos k y ) and 
e Tik = — ti[l + e lk * + e lky + e 2 ( fc *+ fe y)]. Here, we have set 
the lattice constant a — 1. In deriving Eq. (2), we have 
introduced c A ( B ) a ,ija = ^Ek c -4(B)a,k ff e l(S; " :r ' + ' :! ' ! ' i) 
and have taken the canonical transformation 
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where r n k = 



(3) 

e 7\k e T,k) and fl «,k = \{tA.k - 

es.k) + (-l)"yj(eXk ~ £ s,k) 2 + eT,k£T,k- 

Eq. (2) describes analytically four energy bands with 
the indexes < u,v >. In Fig. 2, we plot these bands along 
the path (0,0) — > (tt, 0) — > (tt, 7t) — > (0,0). In our calcu- 
lations, we have used ii = 0.5,^2 = 0.2, t% = — 1.0, t^ — 

0. 02, and fi = -0.622 (half filling) (eV). Obviously, there 
exist two hole Fermi surfaces around (0,0), i.e. a- and 
/3-bands, and two electron Fermi surfaces around (tt,tt), 

1. e. 7- and 5-bands. This is consistent with those ob- 
served from ARPES experiments [7-15]. We note that 
the hole and electron pockets are associated with u — 1 
and while v — and 1 represent the inner and outer 
Fermi surfaces of the hole and electron pockets, respec- 
tively. The parameters t\ , t% and <3 determine the sizes 
of the hole and electron pockets, and ti controls the in- 
tervals between the inner and outer Fermi surfaces. We 
note that /i < —0.622 and /i > —0.622 correspond to hole 
and electron dopings, respectively. When hole (electron) 
doping increases, the hole (electron) Fermi surfaces, i.e. 
a- and /3-bands (7- and (5-bands), become larger while 
the electron (hole) Fermi surfaces i.e. 7- and <5-bands 
(a- and /3- bands) become smaller. The variation of the 
Fermi surfaces with hole or electron doping has been ob- 
served by ARPES experiments [7-15]. When fj, > -0.32 
(i.e. ~ 26.5% electron doping), /3-band disappears. 

In order to investigate superconductivity in iron ar- 
senides, we now introduce the mean field BCS Hamilto- 
nian 
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FIG. 3: (Color online) (a) The differential conductance as 
a function of the bias voltage eV for the pairing symmetry 
A^^ = | Ao(cos k x +cos k y ) with Ao = 5.8 meV and Au.^l 
at optimal electron doping (~ 15%) under temperature 4.2K. 
(b) The corresponding Fermi surfaces. 

where A u „ : k are the superconducting gaps on the energy 
bands < it, v >, depending on the momentum of the long- 
lived quasiparticles ip U v,ka-- We can see from Eqs. (3) and 
(4) that both inter- and intra-band pairings in the orig- 
inal electron operators CAa,ka and CBa,ka- are automati- 
cally included. Here, we assume that the pairing poten- 
tial between electrons is unique, which can avoid many 
superconducting transition temperatures [33]. In other 
words, the energy gaps on all the Fermi surfaces can be 
fitted by a single function of the momentum, i.e. Aio,k = 
An t k = A o,k = A i,k- In Ref. [12], Nakayama et al. 
measured the energy gaps on different Fermi surfaces 
in optimally hole-doped Bao.6Ko.4Fe2As2 (T c ~ 37K) 
by employing ARPES experiments. The order param- 
eter can be fitted as A^k = \^o{cosk x + cosk y ) with 
Ao = 13.5 meV or |A utI .k|. However, in the STM ex- 
periments on optimally electron-doped BaFe1.sCoo.2As2 
(T c ~ 22. 5K) [31,32], only two coherence peaks were ob- 
served at a small gap, i.e. ~ ±5.8 meV. In the following 
we shall calculate the differential conductance for the ex- 
tended s-wave pairing symmetry in the optimal electron 
doping, so that we can compare our theory with the STM 
experiments. 

After diagonalizing the Hamiltonian H = Hq + Hsc 
by the Bogoliubov transformation, we obtain the local 
density of states (LDOS) on the sublattices A or B 

a A A ' B £ 2 
iVvr ^-f iu n - {-lyvluvM 

UVK.V 

(5) 

where v = ±1, < k = < k At k , < k = W T J V l.^ 

= ^{E UV M- f i) 2 + A 2 uv k , and ^ )k „ = \[l + 

(— l) v E Q ,k ~* 1 ]- Obviously, the quasiparticles on the hole 

and electron pockets have different weights A u k to con- 
tribute to the LDOS. 

Usually the STM experiments are performed at low 
temperatures. In order to compare accurately with the 
STM experiments, we must take the effect of temperature 
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FIG. 4: (Color online) The LDOS near an impurity, which 
includes a dominant nonmagnetic potential V s and a small 
magnetic part V m , for the pairing symmetry A uv> ^ = 
^A[}(cosk x + cosky) with Ao = 5.8 meV and lA^kl at opti- 
mal electron doping (~ 15%). 

into account. The differential conductance measured by 
the STM experiments is 

dl f°° 

— oc - J f'(uj - eV)p A ,B( UJ )du, (6) 

where /' is the derivative of the Fermi function and V 
is the bias voltage applied between the STM tip and the 
sample. 

According to the formulas (5) and (6), we can calcu- 
late the differential conductance with different pairing 
symmetries and dopings at low temperatures. In Fig. 
3(a), we present the differential conductance for the ex- 
tended s-wave symmetry with optimal electron doping 
under temperature 4.2K. We have observed that the main 
contribution to dl/dV comes from the hole Fermi sur- 
faces, i.e. a- and /3-bands. Therefore, whether or not 
the nodal points on 7- and (5-bands exist does not change 
qualitatively the features of dl/dV. The main difference 
between the theoretical results and the STM data could 
be due to the fact that either (3- or 5- band of the STM 
sample is much closer to the nodal line as depicted in 
Fig. 3(b). 

In order to detect the sign reversal pairing in the 
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FIG. 5: (Color online) The LDOS near a unitary impurity 



We also investigate the other cases of the impurity 
potential. For the attractive scattering potential, i.e. 
V s < 0, the stronger resonance peak inside gap in the 
LDOS near the impurity site always appears at negative 
energy. With increasing \V S \, the resonance peaks be- 
come higher. When — > — oo, the LDOS is identical 
with that for V s — > +oo, shown in Fig. 5. We note that 
an extra small magnetic potential does not change the 
features of the LDOS. For a dominantly magnetic po- 
tential, the LDOS near the impurity site for A M „ ; k has 



for the pairing symmetry A„„, k = ±A (cos k x + cos k y ) with similar structures with that for |A u „, k |, although the val- 



Ao = 5.8 meV and 
15%). 



| Ami,, k | at optimal electron doping 



Fe As-based superconductors, now we calculate the LDOS 
near an impurity located at the origin of the sublat- 
tice A described by H lmp = V s J2 aa c Aa.ooa c Aa,oacr + 

ooT CAct <ooT c Ua, oo i CAa -. oo l) - Here, V s and V m 
represent the nonmagnetic part and magnetic part of the 
impurity potential, respectively. The total Hamiltonian 
H = Hq + Hsc + Hi m p can be solved by T-matrix ap- 
proach [34]. The analytical expression for the LDOS on 
the sublattices A and B near the impurity has been de- 
rived and will be presented elsewhere [35] . We note that 
the inter-band scattering is only allowed for those bands 
with the same index v. 

In Fig. 4, we plot the LDOS curves for A u „,k and 
|A M1 ,,k| on and near the impurity site with a moderate 
strength of nonmagnetic potential, i.e. V s — 0.25 eV, 
plus a small magnetic potential, i.e. V m = 0.08 eV. Ob- 
viously, for a pure scattering potential (V m = 0), the 
LDOS for A ul ,.k has two impurity resonance peaks at 
±ujq = ±2.8 meV on the impurity site and has a sharp 
resonance peak at ujq near the impurity site. However, 
the LDOS for |A MI)i k| has no such in-gap impurity states. 
Therefore, these impurity resonances can be used to de- 
tect the sign-reversal pairing in the FeAs-based supercon- 
ductors. The origin of these impurity resonances comes 
from the Andreev's bound states due to the inter-band 
quasiparticle scattering with the phase opposite order pa- 
rameters, similar to that of the zero bias resonance peak 
on the Zn impurity in cuprate superconductors. An ad- 
ditional small magnetic potential can strongly suppress 
the impurity peak at — luq and enhance the impurity peak 
at (x>o on the impurity site. Meanwhile, all the resonance 
peaks on different sites slowly move forward to zero en- 
ergy. We note that for the mixing potential, the LDOS 
for |A ut ,.k| is similar to that induced by a magnetic im- 
purity in s-wave superconductors. 

Fig. 5 shows the LDOS for A M1)j k and [A^kl near 
the impurity site with a unitary potential. The LDOS 
for A ul>i k also has two impurity resonance peaks at ±ljq. 
However, the resonance peak at — luq is much stronger 
than that at ojq. For |A mi k|, the LDOS also has no in- 
gap impurity resonance peaks. 



ues of V m or the locations of in-gap resonance peaks are 
different. Therefore, the magnetic impurity seems not to 
be good tool to detect the sign reversal pairing in the 
FeAs-based superconductors. 

In summary, we have built a two-orbital four-band 
tight-binding model describing correctly the character- 
istics of the Fermi surfaces in the FeAs-based supercon- 
ductors. In the framework of mean field theory, we have 
studied the differential conductance and the impurity ef- 
fect for extended s-wave pairing symmetry. It is shown 
that the in-gap impurity resonances induced by nonmag- 
netic scattering potential can be regarded as a signature 
of sign-reversal pairing symmetry in the FeAs-based su- 
perconductors, which could be detected by STM exper- 
iments. These resonance peaks also exhibit in the over- 
doped and underdoped FeAs-based superconductors. 
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